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Abstract 

We calculate the next-to-leading order perturbative corrections to the SVZ sum rules for the coupling f^, the nucleon 
leading twist wave function at the origin. The results are compared to the established loffe sum rules and also to lattice 
QCD simulations. 



1. Introduction 

Hard exclusive reactions have long been recognized as an important tool in the exploration of the nucleon structure 
at different scales. The emergence of quarks and gluons as the adequate degrees of freedom is expected to happen at 
momentum transfers that are accessible in present and planned experiments, in particular on nucleon electromagnetic 
form factors and electroproduction of nucleon resonances see e.g. Uj-^J- 

Quantum chromodynamics (QCD) predicts Isj-lst] that at large momentum transfer the form factors become in- 
creasingly dominated by the contribution of the valence Fock state with small transverse separation between the par- 
tons. In the collinear approximation the wave function can be written in terms of the momentum fraction distribution 
of the three valence quarks 0-i{xi,x2, X3), dubbed distribution amplitude 



[dx\ 

-.<P-i{XuX2,X-i) 



x^'''ul^{x,)[u'^^{x2)dl^(x,) - dl^{x2)ul^{xi)] 10) , 

where the integration measure for the momentum fractions is defined as {dx\ - dx\dx2dxT,6{x\ + X2 + xt, - 1) and 
arrows indicate quark helicities. The integral of the distribution amplitude defines a dimensionful constant 



In^J [dx] 02(xi,X2,X3) , (2) 

which determines the nucleon wave function at the origin. It is a fundamental (scale-dependent) nonperturbative 
constant which plays the central role in QCD description of hard exclusive reactions with protons, and determines 
an overall normalization of the amplitudes. In the academic limit of very large momentum transfers the shape of the 
distribution amplitude 02(xi,X2, X3) is fixed by asymptotic freedom ||7l|81] and e.g. the proton magnetic form factor is 
determined entirely by fff. 

This constant has been estimated several times in the past using QCD sum rules ll9l-[Tlll and more recently also in 
lattice QCD simulations fl2]. Different QCD sum rule calculations are consistent with each other but the lattice result 
appears to be 30% lower This discrepancy calls for a reevaluation of QCD sum rules for including higher-order 
(NLO) contributions, which is the goal of this work. Another motivation to derive the QCD sum rule for /jy t o NLO 



accuracy is that it enters calculations of baryon form factors using the light-cone sum rule approach 0131 11411 which 
are currently being advanced to the NLO as well ifTsll . 

The presentation is organized as follows: In Sec.|2]we introduce two relevant interpolating currents for the nucleon 
and define the corresponding normalization constants. The general properties of their two-point functions as well as 
the diagrams contributing to them are presented in Sec. [3] The results of our calculations are listed in Sec.|4] followed 
by a numerical analysis in order to extract a new value for in Sec. |5] We conclude in Sec. |6l 
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2. Proton interpolating currents 



The leading twist normalization constant can be defined via a nucleon matrix element of a local three-quark 
operator: 

riczix) = \e'''''\{u"\x)Ct.u\x))y,td'{x) 

3 '^^ ' (3) 

-[u'''^{x)Ctd''{x))yit'^'{x)\ . 

Here, a,b,c are color indices, C is the charge conjugation matrix {CyJ^C^^ - -y^i) and z is an arbitrary light-cone 
vector (z^ - 0). The prefactor 2/3 has been chosen to obtain the normalization given below in (ISal i. Note that rjcz is 
constructed such that only isospin-1/2 contributions emerge when acting on a baryon state. In what follows we refer 
to (O as the (isospin improved) Chernyak-Zhitnitsky (CZ) current |9]. 



For comparison we will also consider sum rules for another operator, known as loff'e current 01611 : 



mix) = e"'''{u''^{x)Cy^u''{x))ysy''d--\x) . (4) 

The matrix elements of these operators between vacuum and the proton state define the coupling constants, in standard 
notation: 

(01 77cz(0) \P(q)} = fNiqzW , (5a) 
{0\i]m\P(q))^AimNN . (5b) 

Here A' is a nucleon (proton) spinor and niff ^ 938 MeV the nucleon mass, - m^. In our calculations the three quark 
fields in the currents will be taken to be massless. For an analysis of the properties of baryonic currents containing 
two massless and one massive quark see iITtIi . 

To avoid confusion we note that the choice of the current is determined by the purpose of the calculation. In 
most QCD sum rule calculations of nucleon properties the loffe current is usually adopted as a standard choice since 
it produces stable results. In these applications the value of the coupling Ai is not interesting by itself and usually 
cancels out in final results (e.g. in sum rule ratios). In our case, however, it is the coupling itself which is of interest, 
so that we are bound to use the CZ current. 



3. Two-point functions 

The starting point for our calculation is a two-point function n(q) for a generic baryonic current ijix): 



n(q) = / J d'x e'"- (01 T{;7(x)77(0)) 10) , (6) 

with T denoting a time-ordered product and |0) the vacuum. 

The two-point function can be decomposed into two parts via 

n{q) = ^77i(<?2) + 1772(^2) (7) 

From dimension counting it becomes clear that 77i can only contain contributions from even-dimensional operators, 
whereas all odd-dimensional ones are contained in /Ta- 

We will consider the following contributions to Ui : The unity operator, the two-gluon condensate and four-quark 
condensates. They will be calculated in perturbation theory up to next-to-leading order. For the loffe current the next- 



to-leading order contributions have already been calculated, e.g. in 11181 [19] and the corresponding sum rules have 
been analyzed in |20]. For the Chernyak-Zhitnitsky current however, only the condensates are known at 0(as) ifioll . 
We will try to improve the situation by calculating all contributions at that order and presenting the results below. 

We performed all calculations assuming massless up- and down-quarks and using dimensional regularization in 
the MS-scheme. Diagrammatical representations for all contributions to the two-point function (|6]l of two Chernyak- 
Zhitnitsky currents (O are given in Fig.[Tl 
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(d) (e) (f) (g) 

Figure 1: Feynman diagrams contributing to the CZ two-point function at 0(as). 



4. Sum rules 



In order to extract values for the normalization constants we consider Borel-transformed sum rules 112 111 . Sum 
rules are a standard technique whose details will not be explained here, see f23\ for an introduction to the method. In 
Tab. [U we show detailed results where the common factor {qz)jf. has been removed and the Borel-transformation 



-q-/n=M- 



dq- 



2 



(8) 



has been performed. Each entry corresponds to a diagram shown in Fig. [T| 

All four-quark condensates have been reduced to the square of the two-quark condensate using the factorization 
hypothesis 12111 (with superscript color- and subscript Dirac-indices): 



{qqY 



(9) 



{4Ncy- 

The calculated condensate contributions have been found to be in agreement with 

Using the results from Tab.[T]and the hypothesis of quark-hadron duality, which introduces the effective threshold 
So, we postulate the sum rule for in next-to-leading order including operators up to dimension 6. The full sum rule 
reads: 



2\fyf e--i 



s e «- 











;r\180 





which constitutes the main result of this work. In addition we also present the sum rule obtained from ||2 



2{2nf\Xi\^ ml e- 



1 53 



TT \12 



In 



+ -M'-Ei + -a' 
4 3 



aJ5 1/ So „ 



1 m| 
3M2 



(10) 



(11) 
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Table 1: Contributions to the sum rule for CZ-currents. 
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The standard abbreviations used here are a 



-(Infiqq), b = {2nf {^G% ml = ^ (0.65-0.8) GeV^, 

E„ - \ - e^M^ 2"rQ Tif-p-) ■ The function & is defined as &{x) - For better comparison with ( fTOl i we can 

rewrite the first part of (1111) in an integral form: 



2(27r)Vil'«4 



i 1 r 

2 J 



ds e 



i.^fZi. 

;r\12 



ln4 



4 



ff./S ^ 1/, So ^ „/ 50 \\\ 1 ml 



(12) 



5. Numerical analysis 

We can now proceed to analyze the new sum rule (fTOl i. For comparison we will also consider the next-to-leading 
order sum rule for the loffe current which has been derived in [20]. 

Let us consider the normalization constants obtained from the sum rules both at leading order and at next-to- 
leading order, shown in Fig.|2]and Fig. [3] respectively, as a function of the Borel parameter M^. The numerical values 
of the other parameters have been chosen as follows: The renormalization scale has been fixed at = 1 GeV ^ m^. 
For the value of the quark condensate we used a recent determination (/U = 1 GeV) ll23ll 



{qq} = -(242 ± 15MeV)^ 



(13) 



while a standard value has been used for the gluon condensate: 



I^G^ ) = (0.012 ± 0.006) GeV^ . 



(14) 



The threshold parameter has been varied around ^/sq ^1.5 GeV, which is a standard choice due to being of the same 
size as the mass of the lowest hadronic resonances. The leading order sum rules have been found to be quite stable in 



4 



0.0065 r 



0.0045 



0.0040 




0.0 0.5 1.0 15 2.0 2.5 3.0 3.5 4.0 

M ^ [GeV^] 

(a) Analysis of the sum rule for 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

r 2 r/--.ir2i 



M ' [GeV^] 
(b) Analysis of the sum rule for \Ai\ 



0.175 




0.160 

0.0 05 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

M ^ [GeV^] 

(c) Analysis of the ratio |/Af|/|ii| 

Figure 2: Plots for the coupling constants ff^ and Ai in leading order at sq = (1-3 GeV)^(short dashed), ( 1 .4 GeV)^(long 
dashed), (1.5 GeV)2(solid), (1.6 GeV)2(dotted) and (1.7 GeV)2(Ught dotted). 
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Figure 3: Plots for the coupling constants in next-to-leading order at sq - (1.3 GeV) (short dashed), (1.4 GeV) (long 
dashed), (1.5 GeV)2(soUd), (1.6 GeV)2(dotted) and (1.7 GeV)2(light dotted). 
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Table 2: Values for the normalization constants obtained from the sum rules ( fTOl i and (fTTI) in leading and next-to- 
leading order compared to lattice results from fl?) (at = 1 GeV). 





IM [GeV2] 


1^1 1 [GeV2] 


IM/I^il 


LO 


(4.7 + 0.7) ■ 10-^ 


(2.8 ± 0.6) ■ 10-2 


1.7-10-1 


NLO 


(5.1 +0.8) -10-3 


(3.4 ±0.8) -10-2 


1.5-10-1 


Lat 


(3.2 ±0.2) -10-3 


(3.6 ±0.2) -10-2 


0.9-10-1 



this region, though the next-to-leading order sum rules turn out to be a little less stable, cf. Fig.|2]and Fig. [3] As one 
can see from the plots the next-to-leading order correction is positive for both | and \Ai\. 

Tab. |2] summarizes the various values for the normalization constants. For the sum rules the values have been 
extracted from the Borel window 1 GeV^ < < 3 GeV^. The errors are estimated by varying the nonperturbative 
parameters according to (fT3l l and ( fT4l l. incorporating a 10% uncertainty for □'5(1 GeV) ^ 0.5 and also varying the 
threshold yfsQ by 10% (i.e. adjusting io by 20%). Other sources of error such as NNLO corrections, operators 
of higher dimension, condensate factorization, etc. are not easily estimated and thus not included. For the lattice 



simulation the quoted errors cover the combined statistical and systematic uncertainties given in II12I1 . 

It can be observed that our NLO correction to Ifyl amounts to ~ 10%, while | is increased by as much as ~ 20%. 
It should be noted that the new value for of 5.1 - 10^ GeV^ does not coincide with recent lattice simulations on 
that matter, which put the value at about 3.2 - 10'^ GeV^. 

The discussion above only covered the magnitudes of the normalization constants. Their phases cannot be fixed, 
since the sum rules only depend on the square of the absolute value. However, we can make a statement about the 
relative phase between and A]. By considering a non-diagonal correlator (between rjcz and rjj) and taking a ratio 



of sum rules it can be determined that the ratio /n/Ai is a negative real number Il24ll . A standard choice is to designate 
as positive and Ai as negative. 

6. Conclusion 

The next-to-leading order perturbative correction to the two-point function of leading-twist currents has been 
calculated in this work. Based on this correlator an improved SVZ sum rule was proposed and used to obtain a new 
value for the leading-twist normalization constant /ai. The leading order result is in agreement, within errors, with 
previous determinations ^ [loi [l4l] . the next-to-leading order corrections increase the size of fy by approximately 
10%. 

The correction turned out to be less significant than in the case of the coupling constant Ai which is associated 
with the loffe current. For this current the next-to-leading order corrections had been calculated before and the 
established sum rules have been analyzed in comparison. There it has been observed that the 0(as) contributions raise 
the magnitude of Ai by approximately 20%. 

The discrepancies between the sum rule estimations and lattice simulations could not be resolved. The calculated 



correction causes an increase of fy, whereas the results obtained on the lattice 111211 are lower than the traditional sum 
rule predictions. However, the lattice results have not yet stabilized and suffer from uncertainties. These originate, 
for instance, from chiral extrapolation, which is due to the fact that the pion-mass dependence of is not fully 
understood. 
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